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Abstract 

The differential cross section for the process 77 tt^tt^ has been measured in the kinematic 
range 0.6 GeV < < 4.1 GeV, | cos0*| < 0.8 in energy and pion scattering angle, respectively, in 
the 77 center-of-mass system. The results are based on a 223 fb~^ data sample collected with the 
Belle detector at the KEKB e'^e~ collider. Using the data with W > lA GeV, we obtain results 
on light-quark resonances and charmonia. We also compare the observed angular dependence and 
ratios of cross sections for neutral-pair and charged-pair production to QCD models. Differential 
cross sections are fitted in the energy region, 1.4 GeV < W < 2.2 GeV, with a simple model 
where partial waves consist of resonances such as /2(1525), /2(1950) and /4(2050) and smooth 
backgrounds. In the higher energy region, we observe production of the XcO charmonium state and 
obtain the product of its two-photon decay width and the branching fraction to The energy 

and angular dependences above 3.1 GeV are compatible with those measured in the tt~^tt~ channel, 
and in addition we find that the cross section ratio, (7(7r°7r'^)/(T(7r"'"7r~), is 0.32 ± 0.03 ± 0.05 on 
average in the 3.1-4.1 GeV region. 

PACS numbers: 13.20.Gd, 13.60.Le, 13.66.Bc, 14.40.Cs,14.40.Gx 
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I. INTRODUCTION 



Measurements of exclusive hadronic final states in two-photon collisions provide valu- 
able information concerning physics of light and heavy-quark resonances, perturbative and 
non-perturbative QCD and hadron-production mechanisms. So far, we have measured the 
production cross sections for charged-pion pairs [l|, 0] , charged- and neutral-kaon pairs [21, 0] , 
and proton-antiproton pairs We have also analyzed D-meson-pair production and ob- 
serve a new charmonium state [5|. Recently, we have presented a measurement of neutral- 
pion pair production based on a data sample corresponding to an integrated luminosity of 
95 fb~^ [6|. We have carried out an analysis to extract information on light quark resonances 
from the energy and angular dependences of the differential cross sections (DCS), by fitting 
to the resonance parameters of /o(980), /2(1270) and other hypothetical resonances. 

Here we present a measurement of the DCS, da/d\ cos^*|, for the process 77 tt^tt^ in 
a wide two-photon center-of-mass (cm.) energy {W) range from 0.6 to 4.1 GeV, and in the 
cm. angular range, | cos^*| < 0.8. We use a 223 fb~^ data sample, which is more than twice 
as large as that in our previous analysis jl, 0|. We focus on the range W > lA GeV, where 
the previous data was statistically limited. 

In the intermediate energy range (1.0 GeV < W < 2.4 GeV), the formation of me- 
son resonances decaying to vrvr is the dominant contribution. For ordinary qq mesons in 
isospin conserving decays to vrvr, the only allowed I^J^^ states produced by two photons 
are 0+(even)++, that is, fj=even mesons. Several mesons with these quantum numbers are 
suggested by results of hadron-beam or charmonium decay experiments in the 1.5 - 2.2 GeV 
region, but none of them have been firmly established in two-photon processes, which are 
sensitive to the internal quark structure of the meson. In addition, the tt^tt" channel has 
two advantages in the study of resonances: a smaller contribution from the continuum is 
expected in it than in the tt+tt" channel; and the angular coverage is larger (| cos6'*| < 0.8 
instead of 0.6). 

At higher energies, we can invoke a quark model. In leading order calculations j^, 0| , which 
take into account the spin correlation between quarks, the vr^vr'' cross section is predicted to 
be much smaller than that of vr+vr", and the ratio of 7r°7r° to vr+vr" cross sections is around 
0.03-0.06. However, higher-order or non-perturbative QCD effects can modify this ratio. 
For example, the handbag model, which considers soft hadron exchange, predicts the same 
amplitude for the two processes, and thus the expected ratio is 0.5 [lO|. Analyses of energy 
and angular distributions of the cross sections are essential for determining properties of the 
observed resonances and for testing the validity of QCD models. 

The organization of this article is as follows. In Section [Til a brief description of the Belle 
detector is given. Section UTTl explains the procedure used to obtain differential cross sections. 
Section [IV] is devoted to results on resonances obtained by fitting differential cross sections 
in the range 1.4 GeV < W < 2.2 GeV. Section |V] describes analyses at higher energy. The 
topics included there are the angular dependence as a function of W, the charmonia Xco and 
Xc2 states and the ratio of tt^tt^ to tt+tt" cross sections. Finally, Section IVII concludes this 
report. All of the results presented in this report are preliminary. 



II. EXPERIMENTAL APPARATUS 



We use a 223 fb 

tor [ii 



data sample from the Belle experiment [ll[ at the KEKB accelera- 
The data were recorded at e+e~ cm. energies of 10.58 GeV (179 fb"^), 10.52 GeV 
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(19 fb-i), 10.36 GeV {T{3S) runs, 2.9 fb-^), 10.30 GeV (0.3 fb"!) and 10.86 GeV (T(55) 
runs, 21.7 fb~^). The difference of two-plioton fluxes per e^e~-beam luminosity (=luminosity 
function) in the measured W regions due to the difference of the beam energies is small (max- 
imum ± 4%). We combine the results from the different beam energies. The effect on the 
cross section is less than 0.5%. 

The analysis is carried out in the "zero-tag" mode, where neither the recoil electron nor 
positron are detected. We restrict the virtuality of the incident photons to be small by 
imposing strict transverse-momentum balance with respect to the beam axis for the final- 
state hadronic system. 



A comprehensive description of the Belle detector is given elsewhere [ll[ . We mention here 



only those detector components that are essential for the present measurement. Charged 
tracks are reconstructed from hit information in a silicon vertex detector and a central 
drift chamber (CDC) located in a uniform 1.5 T solenoidal magnetic field. The detector 
solenoid is oriented along the z axis, which points in the direction opposite to that of the 
positron beam. Photon detection and energy measurements are performed with a CsI(Tl) 
electromagnetic calorimeter (ECL). 

For this all-neutral final state, we require that there be no reconstructed tracks coming 
from the vicinity of the nominal collision point. Therefore, the CDC is used for vetoing 
events with charged track(s). Photons from decays of two neutral pions are detected and 
their momentum vectors are measured by the ECL. Signals from the ECL are also used to 
trigger signal events. 



III. DERIVING DIFFERENTIAL CROSS SECTIONS 

The event triggers, data processing, and event selection are the same as those described 
in Ref. j^. We derive the cm. energy W of the two-photon collision from the invariant 
mass of the two neutral pion system. We calculate the cosine of the scattering angle of 7r° 
in the 77 cm. frame, | cos^*| for each event, using the e~^e~ collision axis in the e^e~ cm. 
frame as the reference of the polar angle as an approximation, because we do not know the 
exact 77 collision axis. 



A. Data reduction 

We find that the signal candidates in the low energy region {W < 1.2 GeV) are con- 
siderably contaminated by background events. We study the p^-balance distribution, i.e., 
the event distribution in 1 1] | , to separate the signal and background components. We 
estimate the unbalanced background component for W < 1.2 GeV, and subtract the yield 
in the signal region, in the same manner as in the previous analysis However, above 
1.2 GeV, we cannot quantitatively determine the background contamination because of the 
small background rate and low statistics of the sample, as well as the uncertainty in the 
functional form for the signal shape. 

Using the ratio of yields between the pt-balanced and unbalanced regions, we can estimate 
the backgrounds. In Fig. 1, we plot the W dependences of 

R = Yield(0.15 < | J^p*] < 0.20 GeV/c)/Yield(| < 0.05 GeV/c). (1) 
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FIG. 1: The yield ratio R in the pt-unbalanced bin to the pt-balanced (signal) bin (see text 
for the exact definition) for the experimental data. The solid line shows the signal-component 
obtained from the signal-MC and corrected taking into account the poorer momentum resolution 
in experimental data. 

We integrate over all angles in this figure. The main part of the Ty-dependence of R comes 
from the energy dependence of the momentum resolution. The expected ratio from the pure 
signal component (calculated from the signal Monte Carlo (MC) events and corrected for 
the deviation of the pt resolution discussed in the next subsection) is shown by the solid 
line. The excess of R over the line (Ai?) is expected to correspond to the contribution from 
the prunbalanced background. The excess is relatively small above 1.0 GeV, although some 
fine structure is visible there. In the range 1.2-1.5 GeV, AR is undetectably small; for 1.5- 
3.3 GeV, AR ranges between 0.00 and 0.08; above 3.3 GeV, AR is the range from 0.08 to 
0.2. From the R values, we estimate that the background contamination in the signal region 
is ~ -R/4, which is smaller than 3% for 1.1 - 3.3 GeV and around 3% for 3.6 - 4.1 GeV. We 
subtract 3% for 3.6 - 4.1 GeV, and assign a 3% systematic error from this source for the 1.5 
-4.1 GeV range. 

We estimate the invariant-mass resolution from studies of signal-MC and experimental 
events. We find that an asymmetric Gaussian function with standard deviations of 1.9%H^ 
and 1.3%iy on the lower and higher sides of the peak, respectively, approximates the smear- 
ing reasonably well. Based on this information, unfolding is performed using the singular 
value decomposition (SVD) algorithm [1^ at the yield level [6'], and is applied to obtain the 
corrected W distribution in the 0.9 - 2.4 GeV region, using data with observed W values 
between 0.72 and 2.68 GeV. For lower energies, W < 0.9 GeV, the effect of the migration 
is expected to be small because the invariant- mass resolution is small compared with the 
bin width. For higher energies, W > 2A GeV, where the statistics is relatively low and 
unfolding would enlarge the errors, we rebin events into 100 MeV bins without unfolding. 
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B. Calculation of differential cross section 



We determine the trigger efficiency for signal using the detector and trigger simulators 
applied to the signal MC events. The signal MC data for e~^e~ — > e"*"e~7r°7r° are generated 



using the TREPS code [IJ] for the efficiency determination, isotropically in | cos 6**1 at 58 
fixed W points between 0.5 and 4.5 GeV, The angular distribution at the generator level 
does not play a role in the efficiency determination, because we calculate the efficiencies 
separately in each | cos6'*| bin with a 0.05 width. Samples of 4 x 10^ events are generated at 
each W point. Two sets of different background conditions, which were extracted from the 
beam collision data are embedded in the signal MC data in the detector simulation; these are 
put through the trigger simulator and the event selection program. To minimize statistical 
fluctuations in the MC calculation, we ffi the numerical results of the trigger efficiency to a 
two-dimensional empirical function in {W, \ cos^*|). 

The efficiency calculated from the signal MC events is corrected for a systematic difference 
of the peak widths in the p^-balance distributions found between the experimental data and 
the MC events, which is attributed to a difference in the momentum resolution for vr^'s. The 
correction factor is typically 0.95. 



The DCS for each (W, \ cos^*|) point is derived from the following formula: 

da AY - AB 



d\cose*\ AWA\cose*\ J CdtL^^{W)7] 



(2) 



where Ay and AB are the signal yield and the estimated prunbalanced background in the 
bin, AW and A| cos 6**1 are the bin widths, / Cdt and L^^(W) are the integrated luminosity 
and two-photon luminosity function calculated by TREPS [ij], respectively, and t] is the 
efficiency including the correction described above. The energy bin width AW is 0.02 GeV 
for 0.6 GeV < < 1.8 GeV, 0.04 GeV for 1.8 GeV < W < 2A GeV, and 0.1 GeV for 
2.4 GeV <W < A.l GeV. The width of each angular bin is A| cos^*| = 0.05. 

Figure 2 shows the W dependence of the cross section integrated over |cos0*| < 0.8. 
The results are obtained by simply summing da/d\ cos^^*| • A| cos6*\ over the corresponding 
angular bins. We have removed the bins in the range 3.3 GeV < W < 3.6 GeV, because 
we cannot separate the Xco and Xc2 components and the continuum in a model-independent 
way due to the finite mass resolution and insufficient statistics of the measurement. The 
cross section in this region is discussed in detail in Section |Vl 

We show the angular dependence of the DCS at several W points in Fig. 3. Note that 
the cross sections in the neighboring bins after the unfolding are no longer independent of 
each other in both central values and size of errors. 

We estimate the systematic errors for the cross section in each energy bin arising from 
various sources. The systematic errors arise from uncertainties in trigger efficiency (4-30%), 
77° reconstruction efficiency (6%), p^-balance cut (typically 3%), background subtraction 
(0-40%), luminosity function (4-5%), beam background effect for efficiency (2-4%), other 
efficiency errors (4%) and the unfolding procedure (0-4%), for which we show a range when 
the relative error size depends on W. The total systematic error is obtained by adding 
the uncertainties in quadrature and is 10% in the intermediate W region. The systematic 
error becomes much larger at lower W. At higher W, the systematic error is rather stable, 
typically ~ 11%. 
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FIG. 2: The integrated cross section results in the angular regions | cos0*| < 0.8. Data points in 
bins near 3.5 GeV are not shown because of uncertainty from the XcJ subtraction. 
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FIG. 3: The DCS for five selected W points, 0.97 GeV, 1.27 GeV, 1.95 GeV, 2.45 GeV and 
3.15 GeV. 



IV. ANALYSIS OF RESONANCES IN THE RANGE 1.4 GeV <W <2.2 GeV 

Previously, we have obtained a reasonable fit to a simple model of resonances and smooth 
backgrounds in the energy region 0.8 GeV < W < l.Q GeV from the DCS of 77 tt^tt^ 
with a 95 fb~^ data sample [6|. The clear /o(980) peak and the large contribution from the 
/2(1270) can be fitted with parameters determined from tt+tt" data [l|. 

In this section, we extend the analysis to the higher energy region 1.4 GeV < W < 
2.2 GeV (which has some overlap with the previous study) using a higher statistics sample 
of 223 fb^^. It is well known that deriving reliable results on partial waves is difficult. 
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especially in the higher energy region considered here. The goal of this study is to provide 
some information on partial waves in this energy region and to demonstrate the sensitivity 
(and limitations) of our high-statistics data sample (three orders of magnitude more statistics 
than the past experiments) with a large angular coverage (| cos 6**1 < 0.8). 



A. Parameterization of Partial Wave Amplitudes 

In the energy region W < 3 GeV, J > 4 partial waves (next is J = 6) may be neglected 
so that only S, D and G waves are to be considered. The DCS can be expressed as: 



^(77 - ^^°) = \S + ^0 Y,' + Go + Y,' + G, F/f , (3) 

where Dq and Go {D2 and G2) denotes the helicity (2) components of the D and G wave, 
respectively, and Yj" are the spherical harmonics. 

We derive some information on resonances in the DCS fit by parameterizing partial wave 
amplitudes in terms of resonances and smooth "backgrounds". Once the functional forms 
of the amplitudes are assumed, we can use Eq. ([3]) to fit the DCS. From Fig. 10 in Ref. ^ 
(also from the same plot with higher statistics), it appears that the G waves are non-zero 
at ~ 1.8 GeV and are dominated by the G2 wave. Here we assume (and check the 
necessity of) including the /4(2050) in the G2 wave; we note that the /4(2050)'s two-photon 
coupling has not been measured. Since the G2 wave interferes with the D2 wave, and the D2 
wave should contain resonances besides the /2(1270), we include the known resonances, the 
/2(1525) and /2(1950), which are known to couple to two photons (while the /2(2010) does 
not). There are several more resonances that might couple to 77 and vrvr in this mass region 
as listed in Ref. [3]. Thus we assume that the /2(1950) is just an empirical parameterization 
representing these other resonances; we denote it as the "/2(1950)". 

We parameterize the partial waves as follows: 



s 


= %(y)e^*^- + 55, 




Do 


^^l+.o/^^^^^™^^" 


° + Bdo, 


D2 


Wl + .o/^^^^^^°^^^'^ 


+ ^/^ (1525)6 


Go 


= 0, 




G2 


= (2050)6*'^'' + Bg2 





(4) 

where ^/^(y), Af^^i27o), ^/^(i525), Af2{z), ^"/2(i950)" and %(2050) are the amplitudes of the 
corresponding resonances; Bs, Bdo, -Bd2 and Bg2 are "background" amplitudes for S, Do, 
D2 and G2 waves; ro2 is the helicity-0 fraction of the /2(1270); and 4>sy, 'Pdo, (pd2, (p2P, (p2Z, 
029 and 04 are the phases of resonances relative to background amplitudes. An fo{Y) term 
was needed to obtain a good fit in the previous analysis of the energy region 0.8 GeV < 
W < 1.6 GeV [1]. The f2{Z) is needed here to "explain" the dip-bump structure in the 
range 1.5 GeV < W < 1.7 GeV as seen in Fig. 2. We assume that Go = and that G2 
consists only of the /4(2050) and a smooth "background". 
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TABLE I: Parameters of the /2(1525), /2(1950) and /4(2050) [15]. 



Parameter 


/2(1525) 


/2(1950) 


/4(2050) Unit 


Mass 


1525 ± 5 


1944 ± 12 2025 ± 10 MeV/c^ 


Width 


731^ 


472 ± 18 


225 ± 18 MeV 


BiTTTT) 


0.82 ±0.15 


seen 


17.0 ± 1.5 % 


B{KK) 


88.8 ±3.1 


seen 


0.68l°:?^ % 


B{r]ri) 


10.3 ±3.1 


seen 


21 ± 8 % 


Bill) 


1.11 ±0.14 


seen 


unknown 10"^ 



We parameterize resonances with the formula given in Eq. ([5]). The relativistic Breit- 
Wigner resonance amphtude Ar{W) for a spin-J resonance R of mass ttir is given by 



0^0 



, _ 87r(2J±l)m^ ^T,,T. 
^^^^^ "V W ml-W-- zm«r,o, • 



Energy- dependent widths are used for the parameterization of the /2(1270) and /2(1525) [l| 
The resonance parameters given in Ref. |15|] for the /2(1525), /2(1950) and /4(2050) are 
summarized in Tabled Since all (some) of the individual decay fractions of the "/2(1950)" 
(/4(2050)) are not known, we neglect the W dependence of their partial and total widths. 



Background amplitudes are parameterized as follows. 



Bs 


= a,r{W 




-WQ)+Csr + i{as^{W - 


-Wo) 


' + bs^{W- 


Wo)+Csi), 


Bdo 


= aoiW- 


-Wof + bo{W- 


Wo) + Co, 








Bd2 


= a2{W- 


- Wof + h2{W - 


Wo) + C2, 








Bg2 


= agr{W 




- Wo) + Cgr + tiagiiW - 


-Wo 


)^ + bg^iW- 


Wo) + Cg,) 



(6) 

where Wo = 1.2 GeV and we fix Cgr = Cgi = to reduce the number of parameters; floating 
them does not improve fitting. We assume background amplitudes to be quadratic in W for 
all the waves for both real and imaginary parts. The background Do and D2 amplitudes are 
taken to be real by definition. 



B. Fit results 



In the fit, we fix the values of the parameters for the /o(^), /2(1270) and the phases (pdo 
and (j)d2 to those determined in Ref. [l| and j^. The parameters of the /2(1525) are fixed 



to those in Ref. [15|. If the resonance parameters of the "/2(1950)" and /4(2050) are fixed 
at the values given in PDG 15| as summarized in Table [H then the fit is very poor yielding 
X^/ndf = 4.0. Here we quote the results when they are floated. By fixing the phases (pdo 
and (j)d2, the sign ambiguity for B^o and i?D2 can be resolved. 

Here the unfolded DCS are fitted. One to three thousand sets of randomly generated 
initial parameters are fitted for each study. Two solutions (denoted by sol. A and B) of 
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FIG. 4: Fitted curves (solid line) and DCS {d(T/d\ cos0*| (nb)) for VK-bins indicated. 

reasonably good fit quality {x^/ndf = 1.08) are obtained only for the nominal fit, where an 
f2{Z) component is included together with the "/2(1950)" and /4(2050), whose parameters 
are floated. Fit results are shown in Fig. 4 for the DCS and in Fig. 5 for the total cross 
section. Two solutions are indistinguishable. The main difference between sol. A and B is 
the two photon coupling of the "/2(1950)" and f2{Z); that is larger by a factor of 8 (5) in 
sol. B for the "/2(1950)" {f2{Z)), which may be too large. These two solutions arise from 
constructive and destructive interference between resonances and backgrounds. 

Since the two-photon coupling of the /4(2050) has not been measured, a fit without it 
is also given in Table HTl The fit quality is unacceptable, indicating that the /4(2050) has 
a non-zero two-photon coupling. The existence of the f2{Z), a spin-2 resonance with mass 
near 1500 MeV/c^, is controversial. Thus, fits are also made without the f2{Z). There, 
two fits are carried out fixing and floating the branching fraction to two photons of the 
/2(1525). The fitted parameters are listed in Table HTl The fit without the f2{Z) gives an 
unacceptable x^- When the branching fraction of the /2(1525) to two photons is floated 
with no f2{Z) contribution, the value obtained is ten times larger than the nominal value, 
which is unacceptable; thus some resonance hke the f2{Z) is indeed necessary. 



The spin of the "resonance" of mass ~ 1500 MeV/c^ (denoted here as f2{Z) for a spin-2 
resonance) is not known. Thus we also fit by assuming the spin to be denoted as the fo{Z), 
which contributes to the S-wave. In additions, fits with the /4(2050) only, the "/2(1950)" 
only and no new resonances are performed; the results are summarized in Table lllli The 
f2{Z) hypothesis is favored over that of the fo{Z) with about 3a significance, which is 
calculated from the difference of values. 



A study of systematic errors is not performed because we do not know how to estimate 
uncertainty from model dependence. 
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FIG. 5: Total cross section (| cos0*| < 0.8 (nb)) and fitted curves. 



TABLE II: Fitted parameters (1) 



Parameter 




With /2(Z) 




Without /2(Z) 






All included 


no 


^(/^(1525)^77) 


Unit 




sol. A 


sol. B 


/4(2050) 


fixed 


free 




Mass(/4(2050)) 


1935+1^ 


1902+^2 




1865+}^ 


1866 ± 19 


MeV/c2 


rtot(/4(2050)) 


3691^^ 


35011 




3791^^ 


356ti 


MeV 


r^^(/4(2050)) 






(fixed) 


9iti 


68 ±12 


eV 


mass( '72(1950)") 


1852123 


1906+^ 


2191 ±6 


16381? 


174ll?2 


MeV/c2 


rtot('72(1950)") 




394li 




357l|^ 


370l|i 


MeV 






'^-22 


o.il^i 






eV 


mass(/2(Z)) 


15261^ 


15421^ 


1649 ± 5 






MeV/c2 


rtot(/2(^)) 


121 ±9 


217tl^ 


414 ± 15 






MeV 


r^^S(7r°7r°) 




87 ±27 


5031^^ 


(fixed) 


(fixed) 


eV 


(1525) (77) 




1.11 (fi 


xed) 




11.9 ± 1.3 


10-6 



(ndf) 485.0 (450) 485.2 (450) 619.1 (454) 571.0 (454) 517.5 (453) 



V. ANALYSIS OF THE HIGHER-ENERGY REGION 

In general, we expect that theoretical models based on QCD give reasonable predictions 
even for two-photon production of exclusive final-states such as 77 tt^tt^ in the high energy 
region. However, the models do not give information on what energies can be considered as 
high enough. 



The handbag model [10| predicts that the angular dependence of the DCS for 77 — »• 
TTTT goes as ~ sin-^6'*. This prediction is common to charged and neutral pairs of pions. 
Our measurement for the charged-pion process agrees with this expectation above W > 
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TABLE III: Fitted parameters (2) 



Parameter 


With fo{Z) 


/4(2050) 
only 


'72(1950)" 
only 


None 


Unit 


Mass(/4(2050)) 


187611^^ 


1894 ± 7 




- 


MeV/c^ 


rtot(/4(2050)) 


4931^^ 




- 


- 


MeV 


r^^(/4(2050)) 


21311^ 


3lt« 




(fixed) 


eV 


Mass( '72(1950)") 


1752;^^ 




1630 ± 31 




MeV/c2 


rtot( '72(1950)") 


310^26 




362^^1 




MeV 


r^^i3('72(1950)" ^ vrOvrO) 




(fixed) 




(fixed) 


eV 


Mass(/o(Z)) 


1566l||j 








MeV/c^ 


rtot(/o(^)) 










MeV 




ni+52 


(fixed) 


(fixed) 


(fixed) 


eV 


(ndf) 


496.3 (450) 1938.6 (458) 705.0 (458) 2950.0 (462) 





3.1 GeV [2]. The prediction of the cross section ratio, cT(7r'^7r'^)/a"(7r+7r^) is determined by 
isospin invariance to be 0.5. 

However, predictions based on perturbative QCD at leading order suggests that an an- 
gular distribution for will be different from that for 7r~^n . The primary term of the 
leading order (a short-range emission term) for the charged-pion process is well described by 
a sin~^ 6* dependence jsj , while this term vanishes for the neutral-pion process. The shape of 
the angular distribution from the next term (a long-range interaction term) is unpredictable, 
and in general, perturbative QCD models predict smaller cross sections for this next term 
and hence a small neutral to charged cross-section ratio, cr(7r°7r°)/(T(7r+7r^) ~ 0.03. 

A. Angular dependence 

We compare the angular dependence of the DCS in the range | cos6'*| < 0.8 for W > 
2.4 GeV with the function sin~^^*. We also try a fit with an additional term, to quantify a 
possible deviation from sin~^ 9* behavior. We choose this function because it gives relatively 
good fits empirically in a wide range of W. Thus the fit function is parameterized as: 

da/d\cose*\ =a(sin-^r + 6cos^r). (7) 

We fit using a binned maximum likelihood method and 16 bins in the range | cos6*\ < 0.8. 
We know that the effect of charmonia is large in the region 3.2 GeV < W < 3.6 GeV, but we 
cannot separate it in the angular dependence because we cannot assume here any functional 
shapes for the non-charmonium component. The results of the fit for b are shown in Fig. 6 
as well as the fit to the angular distributions in the four selected W regions, where the 
DCS, the vertical axis of this figure, is normalized to the total cross section a{\ cos^*| < 0.8) 
in each W region, i.e. the area under the curve is 1. The parameter b is close to zero 
above > 3.1 GeV (comparing with b = 9.279, which would give the same contribution 
in cr(|cos6'*| < 0.8) as the sin~^6'* term does) although it becomes nearly constant and 
then systematically negative above the charmonium region. The change in the b parameter, 
which approaches a constant value near zero, occurs at a value close to that observed in 
the charged pion case. 
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FIG. 6: (a) The fits of tlie angular dependence of the normahzed DCS (see text) at four selected 
W points. For the blue curves the coefficient h (see the fit formula in the text) is fixed to 0. The 
magenta curves show the fits with h floating, (b) The energy dependence of the parameter h giving 
the best fits. Here, the charmonium contributions are not subtracted, and the data in the XcO ai^d 
Xc2 charmonium regions are plotted with different colors. 



B. Yields of Xcj charmonia 

The structures seen in the yield distribution for 3.3 GeV <W < 3.6 GeV (Fig. 7) are from 
charmonium production, 77 Xco? Xc2 —>■ Tr'^vr". Similar production of the two charmonium 
states is observed in the vr^vr" -, K+K- and TT^K^ final states Q,^. 

We fit the distribution to contributions of Xco, Xc2 and a smooth continuum compo- 
nent, using the following function: 



Y{W) = WakW-P + e^ViV^,,BW^,„(iy)p + N^JB^U^)? + «(1 - (8) 

in the W region between 2.8 and 4.0 GeV, where BW^^^(Vr) is a Breit-Wigner function for 
the charmonium amplitude, which is proportional to ~ — M^^^ — iM^^jT^^j) and is 

normalized as / |BW^^j(iy)p(iiy = 1. The masses and widths, M and F, of the charmo- 
nium states are fixed to the PDG world averages 15|. The component aW~^ corresponds to 
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TABLE IV: Results of the fits (see text) to obtain the charmonium contributions with and with- 
out interference effects. Errors are statistical only. Logarithmic likelihood (ln£) values are only 
meaningful when comparing two or more fits. 



Interference 




(f) N^c2 


-2ln C/ndf 


Without 


100 ±16 




52.4/56 


With 


103t^° 0.82l°i« 


(l.l±0.3)7r 34± 13 


44.2/54 



the contribution from the continuum component, with a fraction k that interferes with the 
XcO ampHtude with a relative phase angle, 0. It is impossible to determine the interference 
parameters for the Xc2, because of its much smaller intrinsic width compared to the mea- 
surement resolution. We fit the Xc2 yield {N^c2) with a formula where no interference term 
is included, and later we estimate the maximum effects from the interference term when 
deriving the two-photon decay width of Xc2- We use data only in the range | cos^*| < 0.4 
where the charmonium contribution is dominant. We take into account the smearing effect 
due to a finite mass resolution in the fit, using the same function as used for the unfolding. 

A binned maximum likelihood method is applied. We examined two cases with and 
without the interference. Reasonably good fits are obtained for both cases. The fit results 
are summarized in Table IIV[ In the table, C is the likelihood value and ndf is the number 
of degrees of freedom. The normalization N^co in Eq-(8) is proportional to the square of 
the resonance amplitude. The yields from the fits are translated into products of the two- 
photon decay width and the branching fraction, T^^{xcj)B{Xcj tt^tt"), which are listed 
in Table El The systematic errors are taken from the changes in the best fits in the central 
values of yields when the absolute energy scale is varied by ±10 MeV for the W measurement 
and the invariant-mass resolution is varied by ±10% for the corresponding Gaussian widths. 
The changes of the goodness of fit (— 21n£) for these variations are found to be small, less 
than 1.7. 

The XcO is observed with a statistical significance of 7.6a {7.3a) when we take (do not 
take) interference into account. The statistical significance for the Xc2 is 2.6a when we take 
interference of the Xco into account, but it is only 1.3a when we do not take into account 
interference. This is because interference makes the line shape of Xco highly asymmetric 
with a short tail and destructive interference on the high-energy side. The red and blue 
curves in Fig. 7 show the fits for the two cases of with and without Xco interference. 

The results for r^^i3(xcj) in the tc^tt^ final state are compared with previous measure- 
ments of r^^S(xcj) in the vr+vr- decay mode, 15.1±2.1±2.3 eV and 0.76±0.14±0.11 eV for 
XcO and Xc2, respectively [2], (or even KK, referring to SU(3) symmetry [3, ^sl]) decay mode. 
Although the effects of interference were neglected in the tt+tt^ measurements, the results 
are consistent with the ratio expected from isospin invariance, B{xcj 7r°7r°)/i3(xcj — *■ 

TT + TT-) =1:2. 



15 



TABLE V: Products of the two-photon decay width and the branching fraction for the two char- 
monia. Here, T^^B{xcj) means T^^{xcj)S{xcj — > vr'^vr''). The first, second and third errors (if 
exists) are statistical, systematic and from the maximal uncertainties of the relative phase in Xc2 
production. 



Interference 


r^7^(Xco) (eV) r^^Bixc2) (eV) 


Without 
With 


9.7 ±1.5 ±1.0 0.18l[J:}^±0.03 
9.9l|:^± 1.0 0.48 ±0.18 ±0.05 ±0.14 




2.8 



3.2 3.4 3.6 

W(GeV) 



3.8 



FIG. 7: The W distribution of the candidate events with |cos0*| < 0.4 near the charmonium 
region. The red and blue curves show the fits described in the text with and without interference 
with the XcO- 



C. Subtraction of the charmonium contributions 

We subtract the charmonium contributions from nearby bins of the charmonium (Xcj) 
region, 3.2 — 3.6 GeV, in order to obtain a pure DCS from the continuum component. We 
use the fit result with interference obtained in the previous subsection. 

The estimated charmonium yield that includes the contribution from the interference 
term is converted to a DCS contribution in each angular bin of | cos 6**1 < 0.8 by assuminga 
flat distribution for the Xco component and a ~ sin^ 9* distribution for the Xc2 component |3| . 
This assumption is only a model. In reality, we do not know how the angular distribution of 
the interference term behaves; the charmonium amplitudes can interfere with the continuum 
components with different J's of unknown sizes. 

For the W = 3.25 GeV bin, the fit result indicates that there is a non-negligible effect 
from XcO when we assume interference, and thus we make a correction for charmonium 
subtraction. The contribution of the charmonium components in the original DCS is 18% 
at I cos ^* I < 0.6. For W = 3.3 — 3.6 GeV, we apply subtraction for the angular bins 
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0.4 < I cos ^* I < 0.8 after extrapolating the charmonium yield determined in the range 
I cosri < 0.4. 

The DCS thus obtained for the continuum is integrated over the range | cos^^*| < 0.6. We 
convert cr(0.4 < | cos6'*| < 0.8) to cr(| cos6*\ < 0.4) for W = 3.3-3.6 GeV, by assuming that 
the angular dependence of the DCS has a ~ sin~^ 9* dependence. The results are plotted in 
Fig. 8. 



D. W dependence and cross-section ratio 

We fit the DCS integrated over angle, a{ \ cos 6**1 < 0.6), to a power law in the cm. energy, 
W~^, for the energy region 3.1 GeV <W < 4.1 GeV, in which the angular dependence of the 
DCS does not show any large changes. In the fit, we do not use the data in the charmonium 
region {W = 3.3 — 3.6 GeV), where we cannot determine the cross section of the continuum 
component in a model-independent manner. 

The obtained result is n = 6.9 ± 0.6 ± 0.7. The systematic error is dominated by the 
uncertainty of the charmonium contribution in 3.1 GeV < W < 3.3 GeV. This value is 
compatible with the results for the tt+tt" and K^K^ processes |2|], but significantly different 



from the case of K'gK'g [3[ 

The fit for 3.1 GeV < < 4.1 GeV is shown in Fig. 8(a). We also show the ratio of 
-K^Ti^ to TT+TT^ cross scctious in Fig. 8(b). The two processes have a similar W^"" dependence 
for 3.1 GeV < < 4.1 GeV, while their ratio is almost constant in this energy region. The 
average of the ratio in this energy region is < Ratio >= 0.32 ±0.03 ±0.05, where the data in 
the 3.3 - 3.6 GeV region is not used when calculating this average. This ratio is significantly 
larger than the prediction of the leading order QCD calculations 0, [9| and is slightly smaller 



than the value of 0.5, which appears in Ref. [lOj based on isospin invariance. 



VI. SUMMARY AND CONCLUSION 

We have measured the process 77 — >■ vr^vr^ based on data from e^e~ collisions correspond- 
ing to an integrated luminosity of 223 fb~^ with the Belle detector at the KEKB accelerator. 
We derive results for the differential cross sections in the center-of-mass energy and polar 
angle ranges, 0.6 GeV < < 4.1 GeV and | cosO*\ < 0.8. 

A simple model is employed to obtain some information on resonances up to 2.2 GeV/c^. 
Differential cross sections are fitted in the energy region, 1.4 GeV <W < 2.2 GeV, with a 
model where partial waves consist of resonances and smooth backgrounds. The G wave is 
seen to be important for > 1.8 GeV, where the G2 wave appears to dominate over the Go 
wave 0. Thus the ^(2050) is included in the G2 wave while Go is set to zero. Parameters 
of the /2(1270) and other resonances are fixed at PDG flE'] values and those determined in 
the previous analysis [l|, [6|. Two solutions with reasonably good fit quality are obtained 
when an additional spin-2 resonance is introduced with a mass near 1.5 GeV/c^ (denoted 
as f2{Z)) in addition to the /2(1525), /2(1950) and ^(2050). Without the ^(2050), the fit 
quality is unacceptable. 

We observe production of the charmonium Xco state and obtain the product of its two- 
photon decay width and the branching fraction to tt^tt". The energy and angular dependences 
above 3.1 GeV are compatible with those measured in the tt+tt^ channel, and we obtain the 
cross section ratio, (T(7r°7r'^)/cr(7r+7r'^), to be 0.32 ± 0.03 ± 0.05 on average in the 3.1-4.1 GeV 
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FIG. 8: (a) The cross sections for the 77 — > vr'^vr'^ (blue diamonds) and 77 — > 7r'^7r~ (violet 
triangles, [2]) for | cos0*| < 0.6. The curve is the fit to the cross section for 77 — > vr'^vr'^ with a 
~ W~'^ functional shape, (b) Ratio of the cross section of the vr^^vr'^ process to the 7r"'"7r~ process. 
The error bars are statistical only. The red line is the average for 3.1 GeV < W < 4.1 GeV. 
The horizontal line (0.5) is an expectation from isospin invariance for a pure / = component. 
In (a) and (b), the estimated charmonium contributions are subtracted in both vr+vr^ and 
measurements. The results in the W region 3.3 GeV - 3.6 GeV (plotted with lighter colors) are 
not used for the fits. 



region. This ratio is significantly larger than the prediction of the leading order QCD 
calculation. 
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